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ABSTRACT 

Gas can be violently stripped from their galaxy disks in rich clusters, and be dispersed over lOOkpc-scale tails or plumes. Young stars 
have been observed in these tails, suggesting they are formed in situ. This will contribute to the intracluster light, in addition to tidal 
stripping of old stars. We want to quantify the efficiency of intracluster star formation. We present CO(l-O) and CO(2-l) observations, 
made with the lRAM-30m telescope, towards the ram-pressure stripped tail northeast of NGC4388 in Virgo. Hll regions found all 
along the tails, together with dust patches have been targeted. We detect molecular gas in 4 positions along the tail, with masses 
between 7x10^ to 2x10® Mq. Given the large distance from the NGC 4388 galaxy, the molecular clouds must have formed in situ, 
from the HI gas plume. We compute the relation between surface densities of star formation and molecular gas in these regions, and 
find that the star formation has very low efficiency. The corresponding depletion time of the molecular gas can be up to 500 Gyr and 
more. Since this value exceeds a by far Hubble time, this gas will not be converted into stars, and will stay in a gaseous phase to join 
the intracluster medium. 

Key words. Galaxies: evolution — Galaxies: clusters: Individual: Virgo — Galaxies: clusters: intracluster medium — Galaxies: 
interactions — Galaxies: ISM 


1. Introduction 

In overdense cluster environments, galaxies are significantly 
transformed, through different tidal interactions, like the ones 
due to other galaxies, the cluster as a whole (e.g. |Merritt|1984| 
Tonnesen et al.|20(j7l l, and the ones with the intra-cluster medium 
(ICM), which strips them from their gas content. This ram- 


pressure stripping (RPS) process has been described by Gunn & 
Gott ( 1972|l and simulated by many groups (|Quilis et al.||2000 


Vollmer et al.||20()T| |Roediger & HensTer||2605f Jachym et al.| 
2007|l. Evidence of stripping has been observed in many cases 
(|Kenney et al.|2004t[Chung et al.|20()7||Sun et al.|2007||Vollmer| 

et al.||2008|l. RPS and/or tidal interactions can disperse the in- 


et al. 2008| l. RPS and/or tidal interactions can disperse the in 

terstellar gas (ISM) of galaxies at large distance, up to lOOkpc 
scales, as shown by the spectacular tail of ionized gas in Virgo 
( [Kenney et al.|20()8] l. 

What is the fate of the stripped gas? According to the time- 
scale of the ejection, the relative velocity of the ICM-ISM in¬ 
teraction, and the environment, it could be first seen as neutral 
atomic gas (Chung et al. [2009 Scott et al. 2012 Serra et al. 
2013jl, then ionized gas detected in Hg (jGavazzi et al.j 200T[ 


Cortese et al.||2007| [Yagi et M^|2007[ [Zhang et al. 


2013 


and 


is finally heated to X-ray gas temperatures (e.g. Machacek et al. 
[2005[[Sun et al.[2010[ ). In rarer cases, it can be seen as dense and 
cold molecular gas, detected as carbon monoxyde (CO) emission 


(Vollmer et al. 2005 1 Dasyra et al.[2012 [Jachym et al. 20141. The 
presence of these dense molecular clumps might appear surpris¬ 
ing, since the RPS should not be able to drag them out of their 
galaxy disks (Nulsen[[l982| Kenney & Young[[1989|l. However 
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they could reform quickly enough in the tail. The survival of 
these clouds in the hostile ICM e nvironment, with temperature 


10' K and destructive X-rays (e.g., Machacek et al. 2004 


Eabian 


et al.|2006 Tamura et al. 1200911 ^ puzzle, unless they are self- 


shielded (e.g. Dasyra et al.|2012| Jachym et al.|2014j l. The pres¬ 
ence of cold molecular gas is also observed in rich galaxy clus¬ 
ters, with cool cores. Here also a multi-phase gas has been de¬ 
tected, in CO, Ha, X-rays and also the strongest atomic cooling 
lines ( [Edge et al.|2010[ ). Ionized gas, together with warm atomic 
and molecular gas and cold molecular gas clouds coexist in spa¬ 
tially resolved filaments around the brightest cluster galaxy, such 
as in the spectacular prototype Perseus A ([Conselice et al.|2001[ 
[Salome et al.|2006[[^TTl[Lim et al.|20T2|. 

The survival of molecular clouds was also observed by 
[Braine et al.[ ( [2000] l in several tidal tails, and in particular in the 
interacting system Arp 105 (dubbed the Guitar), embedded in 


the X-ray emitting medium of the Abell 1185 cluster (Mahdavi 
[et al.|199^. Again, the forma tion in situ of the molecular clouds 
is favored ( [Braine et al. 2000| l. In the Stephan’s Quintet compact 
group, where X-ray gas and star formation have been observed in 
between galaxies ( O’Sullivan et al.|2009[ l, the shock has been so 
violent (1000 km/s) that H 2 molecules are formed and provide 


the best cooling agent, through mid-infrared radiation (Cluver 
[et al.|201^ . In this shock, multi-phases of gas coexist, from cold 
dense molecular gas to X-ray gas. 

Does this gas form stars? In usual conditions, inside galaxy 
disks, the star formation is observed to depend on the amount of 
molecular gas present (e.g. Bigiel et al. [2008[[Leroy et al.[2013| l. 
A Schmidt-Kennicutt (S-K) relation is observed, roughly linear, 
between the surface densities of star formation and molecular 
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gas, leading to a depletion time-scale (xdep = 2igas/2sFR) of 2 
Gyr. But this relation does not apply in particular regions or cir¬ 
cumstances, such as galaxy centers (Casasola et al.|20T5 1, outer 


parts of galaxies and extended UV disks (|Dessauges-Zavadsky 


et al. 

|2014 

2008 

1. Littl 


e is known on star formation in gas clouds stripped 


constraints on this process, concluding to a very low star for¬ 
mation efficiency, lower by an order of magnitude than what is 
usual in galaxy disks, and even lower than outer parts of galaxies 
or in low surface brightness galaxies. It is interesting to better 
constrain this efficiency, given the large amount of intracluster 
light (ICL) observed today (e.g. [Feldmeier et al.|[2002| |Mihos| 
et al.|2003 1. These stars could come from tidal stripping of old 


stars formed in galaxy disks, or also a large fraction could have 
formed in situ, from ram-pressure stripped gas. More intraclus¬ 


ter star formation could have formed in the past (DeMaio et al. 
2015|l. The origin of the ICL could bring insight on the relative 


role of galaxy interactions during the cluster formation, or clus¬ 
ter processing after relaxation. 


1.1. The tail northeast of NGC4388 


One of the environments to probe the survival of molecular gas 
and the efficiency of star formation under extreme ram-pressure 
conditions is the RPS tail north of NGC 4388 in the Virgo clus¬ 
ter south of M86, where X-ray gas has been mapped ( Iwasawa] 
jet al.|2003] l and young stars have been found ( |Yagi et al.|2013| l. 
It is located at about 400 kpc in projection from the cluster cen¬ 
ter M87. NGC 4388 is moving at a relative velocity redshifted 
by 1500km/s with respect to M87, and more than 2800 km/s 
with respect to the M86 group. This strong velocity may explain 
the violent RPS, the high HI deficiency of NGC 4388 ( |Cayatte| 
et al. 1 1990 1 and the large (~ 35kpc) emission-line region found 
by Yoshida et al. (2002 1 , northeast of the galaxy. The ionized gas 
has a mass of 10^ Mq, and is partly excited by the ionizing radi¬ 
ation of the Seyfert 2 nucleus in NGC 4388. The RPS plume is 
even more extended in HI ( Oosterloo & van Gorko m|2005|l, up 


to 110 kpc, with a mass of 3.4 x 10“ Mq. Gu et al. (2013 i have 


found neutral gas in absorption in X-ray, with column densities 
2-3x10^** cm“^, revealing that the RPS tail is in front of M86. 
The high ratio between hot and cold gas in the clouds means 
that significant evaporation has proceeded. Yagi et al. p013| l find 
star-forming regions in the plume at 35 and 66 kpc from NGC 
4388, with solar metallicity and age 6 Myrs. Since these stars are 
younger than the RPS event, they must have formed in situ. 

In the present paper we present CO detections in the ram- 
pressure stripped gas northeast of NGC 4388. In a previous pa¬ 
per, we have already found molecular gas in a ionized gas tail 
south of M86 ( |Dasyra et al.||201^ , and discussed its survival 
conditions. We here study the link between new stars formed and 
molecular gas, to derive the star formation efficiency. In the RPS 
plume, a significant fraction of the Ha emission could originate 


from the ionized gas in the outer layers of molecular clouds ( Per- 
land et al.|2009 1. This makes the Ha lumps good tracers of star 


formation in an RPS tail, to probe the efficiency of the process of 
formation of intracluster stars. Section 2 presents the IRAM-30m 
observations. Section 3 the results obtained, which are discussed 
in Section 4. 


In the following, we assume a distance of 17.5 Mpc to the 
Virgo cluster (Mei et al. 2007). 
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2. Observations and data reduction 


CO observations along the HI plume 
[2005] l connecting NGC4388 and M86 
30-mt telescope at Pico Veleta, Spain, in two separate runs. The 
first run was part of the project 195-13, with 28 hours of ob¬ 
servation, and took place between the 5th and 8th of December 
2013, with excellent weather conditions (t < 0.1 and a pwv be¬ 
tween 0.1 and 3mm). The second run was project 075-14, with 
47 hours of observations between June 25th-30th 2014, and had 
poor to average weather conditions (r between 0.2 and 0.6 and a 
pwv between 3 and 10 mm). 

All observations were done with the EMIR receiver in the 
E0/E2 configuration, allowing us to observe simultaneously 
CO(l-O) and CO(2-l) at 115.271 and 230.538 GHz respectively. 
The telescope half-power beam widths at these frequencies are 
22" and 11" respectively. The observing strategy consisted in 
single ON-hOEE pointings per each target, with wobbler switch¬ 
ing. 


(Oosterloo & van Gorkom 
were done with the IRAM 


Targets along the HI plume were selected for having a match 
of HI (using the Nhi map from | Oosterloo & van Gorkomj: 


Ha (with data from [Kenney et ar[|2008| and | Yagi et al^ 2013i 
and 25Qp emission (HERSCHEL SPIRE data from [Davies et al.| 
[2012[ l. See Eigure[2 With this criteria 6 targets where selected 
for the first run 195-13, and are listed in Table [1] (first 6 rows). 
As a result of this run, only two sources showed CO detection: 
Source-1 and HaR-2. Since HaR-2 is of particular interest for 
being so far away from both galaxies and for having a strong 
Ha detection ( Yagi et al.|2013| l, it was chosen as a central target 
around which other 5 extra targets were selected for the second 
run 075-14 (second half of Table[2l, following the path of an HI 
peak (Eigure[T]top right box) 


Table 1. Targets and Observations 


Source 

RA(J2000) 

DEC(J2000) 

hel. vel. 

ON-hOEE 


I'/ir.m.i j 

^eg 

[km/s] 

[hrs] 

Source-1 

12:26:01.3 

12:42:30.1 

h-2500 

4.2 

Source-2 

12:26:04.9 

12:45:16.7 

h-2500 

2.9 

Source-3 

12:26:25.4 

12:50:53.6 

h-2200 

3.1 

Source-5 

12:26:24.5 

12:58:14.6 

+2000 

3.9 

HaR-2 

12:26:13.7 

12:51:36.9 

+2230 

7.2 

HaR-3-4 

12:26:13.5 

12:43:03.7 

+2500 

3.8 

HaR-2-1 

12:26:16.9 

12:53:55.6 

+2230 

4.5 

HaR-2-2 

12:26:16.9 

12:53:19.3 

+2230 

5.6 

HaR-2-3 

12:26:16.9 

12:52:43.1 

+2230 

6.2 

HaR-2-4 

12:26:15.5 

12:52:01.1 

+2230 

8.0 

HaR-1 

12:26:14.6 

12:51:46.4 

+2230 

3.8 

NGC4388 

12:25:46.6 

12:39:44.0 

+2550 

1.2 

Notes. Sources 

in the first 6 

rows correspond to 

run 195- 

•13, and last 


6 rows to run 075-14 (although source HaR-2 was observed in both 
runs). NGC4388 was observed for calibration purposes. Heliocentric 
velocities are referencial, taken from [Oosterloo & van Gorkom[ j2003) 
along the HI plume. 

Concerning the spectral resolution of our data, during the 
observations both ETS and WILMA backends were used, simul¬ 
taneously. The ETS backend has a spectral resolution of 195 kHz 
and a bandwidth of 32 GHz including both polarizations. At 115 
GHz these values correspond to 0.5 and 83200 km/s, and at 230 
GHz to 0.25 and 41600 km/s. As for the WILMA backend, we 
obtained a spectral resolution of 2 MHz and a bandwidth of 16 
GHz. At 115 GHz this translates to 5.2 and 41600 km/s, and to 
2.6 and 20800 km/s at 230 GHz. 























































































C. Verdugo et al.: CO in the Virgo Cluster 



RA (J2000) 


Fig. 1. Targets observed in the Virgo Cluster. Orange: HI contours levels at 1, 5, 10 and 50 xlO'^cm^^ fro m [Oosterloo & van Gorkom 1 2005| l, 
showing the plume of atomic gas stripped from NGC4388. Purple: Ha contours levels at 5, 11 and 50 jsec ( [Kenney et al.|2008| >. Cyan: 250/jm 
contour levels at 0.01 and 0.1 Jy/beam ( [Davies et al.|2012y In zoomed regions: targets selected close to HaR-2 for run 075-14 (red circles). Circles 
enclosing targets are 22" width, as the CO(l-O) HPBW. 


The data were reduced using the CLASS software from the 
GILDAS package. First, a careful inspection of all scans was 
done, to remove the bad ones. The approved scans of the same 
source, CO line and backend were averaged with a normal time 
weighting, to obtain one spectrum. Then, each spectrum was in¬ 
spected individually, and in both of its polarizations, to identify a 
possible CO emission line. If a detection was found in the spec¬ 
tra of both backends, the best spectrum was chosen (in terms of 
spectral resolution and S/N) as the final one. The selected spec¬ 
tra with CO emission are presented in Figures and which 
contain both polarizations, horizontal and vertical, combined. 

Baselines were subtracted with polynomials of order 0 and 
1 depending on the source, and antenna temperatures were cor¬ 
rected by the telescope beam and forward efficiencie^to obtain 
main beam temperatures. Spectra were smoothed with the han- 


* http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 


ning method to degrade the velocity resolution until obtaining a 
value no greater than 1/3 of the FWHM line. 

Finally, a simple gaussian line was fitted to the line candi¬ 
date. The CLASS fit return the velocity position of the line, its 
FWHM, the peak temperature and the integrated line intensity. 

Such spectra and fitting results for sources with CO detection 
are presented in Figures l^(from the first run) and[^(from second 
run) and then in Table For the rest of the sources, with no 
visible CO detections, 3cr upper limits for Ico where calculated 
from their rms values, assuming a Av = 30 km/s . Such limits are 
presented in Table 

3. Results 

After selecting the final spectra for every source with a visi¬ 
ble detection and fitting their gaussian profiles, line parameters 
are calculated and presented in Table These CO detections 
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Velocity (km/s) Velocity (km/s) 

Fig. 2. Top: Spectra taken at the center of NGC4388, observed for calibration purposes. Bottom: Spectra from observing run 195-13 that showed 
CO emission. Spectra information and line parameters are presented in Table The temperature scale corresponds to main beam temperature in 
mK. 


present velocity centroids in the range of ~2200-2500 km/s, and 
are quite consistent with the HI velocities from Table [T] taken 
from |Oosterloo & van Gorkom||2005 (their Fig. 2). H 2 masses 
were derived from the CO(l-O) line intensity, using a Galactic 
CO conversion factor of 2.0 x 10^°[cm“^(K km/s)“'] and a cor¬ 
rection factor of 1.36 to account for heavy elements: 


molecular, need to be taken into account when estimating a total 
amount of gas to be converted into stars. 

Molecular gas can be directly estimated from the CO(l-O) 
line intensity, obtained from our observations. If we take Eq.[2 
and we divide it by the source area (ttR^), we obtain the H 2 sur¬ 
face density: 


MhJMo] = 4.47rR2[pc] Ico(i- 0 )[K km/s] 


(1) ShJMoPC- 2] = 4.4Ico(i-0)[Kkm s-*] 


( 2 ) 


where the source’s radius R corresponds to the CO(l-O) 
beamsize radius at the distance of the source (17.5 Mpc from 
|Mei et al.|2007l l. 

3.1. Star formation efficiency 

To estimate how fast is the gas being transformed into stars, we 
compare the SFR surface density versus the gas surface den¬ 
sity in a S-K relation, to understand the efficiency of this star 
formation process. Since these are low gas density regions, we 
can expect a large amount of gas in atomic phase, greater than 
in molecular phase. Therefore both components, atomic and 


These values are listed in column 2 of Table|4]for the sources 
with CO(l-O) detections, including upper limits for the sources 
with no CO detection using Table 

For the atomic gas component, we estimated the amount of 
HI from the HI column density map of the NGC4388 plume from 
lOosterloo & van Gorkom| ( |2005| l. The atomic gas mass is derived 
from the integrated amount of Nh inside the source solid angle: 




Mhi = yumn I NndA = I NndQ 


/■ 


(3) 
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C. Verdugo et al.: CO in the Virgo Cluster 




2000 2200 2400 2600 2800 2000 2200 2400 2600 2800 

Velocity (km/s) Velocity (km/s) 

Fig. 3. Final CO spectra for HaR-2, HaR-1 and HaR-2-4 from run 075-14. Only FIaR-2 in CO(2-l) has data combined from both runs. All spectra 
have both polarizations combined. Spectra information and line parameters are presented in Table The temperature scale corresponds to main 
beam temperature in mK. 


Aperture photometry was done in the Nh map of Oosterloo 
|& van Gorkom| ( 2005| ) to obtain the integrated column densities 
for our sources. We used 22" apertures to be consistent with our 
CO(l-O) observations . Since these apertures are smaller than 
the spatial resolution of the HI map (18 x 95.1 arcsec), the pho¬ 
tometry results are equivalent to the pixel value of the HI map at 
the position of our CO targets. These values are listed in column 
3 of Table|^ Then, by dividing Eq.|^by the CO(l-O) beam solid 
angle Q = n^, we obtain the HI surface densities Shi listed in 
column 4 of Table |4] 

Finally, we estimate the SFR surface density Ssfr directly 


from the H„ emission in these regions. From Kennicutt & Fvans 

( |20T2t : 


log SFR[Moyr-'] = log Fnjerg s-'] - 41.27 


(4) 


which, divided by the CO(l-O) beam solid angle O, gives 
the SFR surface density Ssfr- Ho, luminosities for HaR-1, HaR- 
2 and HaR-3-4 were obtained from Yagi et al. p013| l, and are 
listed in Table along with the corresponding Ssfr- For the re¬ 
maining sources we used the Hq, map from|Kenney et al.](|2008|l. 


in diameter, similar to the 
Yagi et ah] ( |2013| l. For the 


using an aperture photometry of ~1' 
seeing of the Hq, observations from 
sources with no visible detection in this map, an upper limit was 
calculated form the noise level of Kenney’s map, estimated in 
4.8 X 10^’ 


erg s 'em ^sr 


The SFR surface densities are plotted as a function of the gas 
surface densities to construct a S-K relation in Figure]^ using the 
values from Table We have plotted separately, the atomic and 
molecular gas component of Egas, along with the combination 
of both, in red, black and blue markers respectively. Arrows de¬ 
note the upper limits values from Table|^ In this Figure, adapted 
from [Jachym et al.| ( |2014| l, we can compare our sources with 
theirs in the Norma Cluster, as well as with the sample of spi¬ 
ral galaxies from |Bigiel et al. ( 2008) 1 (in coloured contours), and 
the sample of 30 nearby galaxies from |Bigiel et al.| ( |201 l| l, plot¬ 
ted as the running medians of Esfr as a function of Ehj, and 
with a typical depletion time of ~2.3 Gyr. Additionally, shaded 
regions have been included to represent regions from the out¬ 
skirts of XUV disk galaxies. NGC4625 and NGC6946 data was 
taken from Watson et al. P014| l (priv. communication), includ¬ 
ing IRAM-30m CO observations, and H„ luminosities measured 
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Table 2. CO detections and their line parameters. 


Source 

line 

backend 

int. time 
[min] 

Vo 

[km/s] 

FWHM 

[km/s] 

resol. 

[km/s] 

Tmb 

[mK] 

rms 

[mK] 

Ico 

[mK km/s] 

Mh, 

[10®h4] 

HaR-2 

CO(2-l) 

FTS 

268 

2224+2 

18+3 

4 

9.8 

2.0 

186+29 


Source-1(1) 

CO(2-l) 

WILMA 

232 

2579+2 

25+5 

5 

6.4 

1.7 

173+29 


Source-1(2) 

CO(2-l) 

WILMA 

232 

2732+2 

16+5 

5 

5.6 

1.7 

99+25 


HaR-2 

CO(l-O) 

WILMA 

318 

2218+4 

26+7 

10 

3.2 

1.5 

88+25 

1.1 


CO(2-l) 

WILMA 

586 

2226+2 

24+6 

10 

6.3 

1.6 

163+29 


HaR-1 

CO(l-O) 

WILMA 

229 

2212+4 

35+13 

10 

5.2 

2.3 

199+54 

2.4 

HaR-2-4 

CO(l-O) 

WILMA 

665 

2404+2 

12+3 

10 

4.0 

1.4 

54+16 

0.7 

NGC4388 

CO(l-O) 

WILMA 

87 

2505+2 

289+5 

5 

54.0 

4.8 

17000+276 

205 


CO(2-l) 

WILMA 

82 

2493+4 

275+7 

5 

67 

8.9 

20000+490 



Notes. First three rows correspond to observing run 195-13. Both components shown in FigureISfor Source-1 are listed. H 2 masses were calculated 
as Mh, [MQ]=4.4;rR^[pc] Ico[K km/s], which uses a Galactic CO conversion factor of 2.0x10^ [cm“^(K km/s)“*] and a correction factor of 1.36 
to account for heavy elements. A distance to the Virgo Cluster of 17.5 Mpc ( |Mei et al.|2007^ was used to calculate the CO(l-O) beam radius as 
source radius. FtaR-2 includes data taken in the first run only for CO(l-O), and data cornbined from both runs for CO(2-l). 

Table 3. CO(l-O) upper limits at 3cr for sources with no detection. 


Source 

int. time 
[min] 

rms 

mK 

T ^ 

ico 

[K km/s] 

Mh, 

[lO^Mo] 

£h2 

[Mopc-2] 

Source-2 

173 

1.2 

<0.11 

<1.21 

<0.48 

Source-3 

186 

1.0 

<0.09 

<0.99 

<0.40 

Source-5 

231 

1.2 

<0.11 

<1.21 

<0.48 

HaR-3-4 

229 

1.0 

<0.09 

<0.99 

<0.40 

HaR-2-1 

269 

1.1 

<0.10 

<1.10 

< 0.44 

HaR-2-2 

334 

1.1 

<0.10 

<1.10 

<0.44 

HaR-2-3 

374 

0.6 

<0.05 

<0.55 

<0.22 


Notes. We assumed a Av of 30 km/s. 


Table 4. Schmidt-Kennicutt relation values. 


Source 

[Mopc-2] 

jNnidQ" 

[lO'^crn-^] 

[Mopc-2] 

log(LHj'’ 

[ergs^'] 

log(ZsFR) 

[Moyr^'kpc^^] 

"^dep(H2) 

[yr] 

Source-1 

0.93 

0.33 

0.40 

<35.51 

<-6.19 

>1.4x10'^ 

Source-2 

<0.48 

0.41 

0.50 

<35.51 

<-6.19 

— 

Source-3 

<0.40 

0.56 

0.68 

<35.51 

<-6.19 

— 

Source-5 

<0.48 

0.60 

0.69 

37.51 

-4.20 

<7.6x10® 

HaR-2 

0.39 

2.12 

2.5 

37.10 

-4.61 

1.6x10'® 

HaR-3-4 

<0.40 

— 

— 

37.75 

-3.96 

<3.6x10® 

HaR-2-1 

<0.44 

5.56 

6.77 

<35.51 

<-6.19 

— 

HaR-2-2 

<0.44 

5.13 

6.24 

<35.51 

<-6.19 

— 

HaR-2-3 

<0.22 

3.67 

4.46 

<35.51 

<-6.19 

— 

HaR-2-4 

0.24 

3.08 

3.74 

<35.51 

<-6.19 

>3.7x10" 

HaR-1 

0.88 

2.05 

2.49 

35.89 

-5.82 

5.8x10" 


Notes. Surface densities consider a solid angle £2=2.74 kpc^, equivalent to the CO(l-O) beamsize at the distance of the source (17.5 Mpc fom Mei 


from 

Oosterloo & van Gorkom 1 2005k. 

® from 

Yagi et al. 

12013k and 

Kenney et al. 

2008 


within a 6" aperture. M63 (NGC5055) data corresponds to the 
bright UV region located at 1.35r25 in |Dessauges-Zavadsky et ah] 
( |2()14| l, with IRAM 30-m CO data, and a SFR measured from the 
FUV and 24/im emission. 

“Isochrones” of constant star formation efficiencies are also 
shown to indicate the depletion times of 10*, 10® and 10'® years 
to consume all the gas, including an additional red isochrone 
to mark the age of the universe as one Hubble time (i.e Tdep = 
13.8Gyr). 
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Contrary to the photometry done in the HI data, the H„ data 
was not measured in a 22" diameter aperture, as the CO(l-O) 
FWHM, but in an aperture of 1" in diameter, similar to the seeing 
of those observations. Since we are averaging this Hq, emission in 
a 22" aperture to calculate the Ssfr. we could be diluting the real 
surface density of the gas being converted into stars. To correct 
for this beam dilution, as a representation we have shifted one 
of the points in Figure to a fictitious corresponding to a 
source’s solid angle of 1" in diameter. This correction translates 




































C. Verdugo et al.: CO in the Virgo Cluster 


O 

Q_ 


>s 

o 

a: 

L. 

00 

w 

cr> 

O 



Fig. 4. S-K relation for sources in Table with filled circles and arrows for values and upper limits respectively. Red markers consider only H 2 
for Egas, black ones only HI, and blue ones the sum of both. Figure adapted from|Jachym et al. j20T4|l, where their sources ES0137-001, 001-A, 
001-B and 001-C have been plotted in a similar way than ours: red circles for H 2 gas, black left arrows for HI upper limits, and blue "error bars" 
to account for both gas components. Colored contours account for the spiral galaxies from |Bigiel et al.| ( |2008) (green, orange, red and purple for 
1, 2, 5 and 10 sampling points per 0.05 dex respectively). Black markers with error bars correspond to the running medians in Esfr as a function 
of (Thj of 30 nearby galaxies f rom|Bigiel et al.| ( |201 l[ l. Shaded ovals represent the data from the outer parts of XUV disk galaxies: NGC4625 and 
NGC6946 from Watson et al.H20l4| (priv. communication) , and M63 (NGC5055) from |Dessauges-Zavadsky et al.|pOT4) , taking only H 2 into 
account in all of them. The dashed vertical line shows the 9 Mopc^^ threshold at which the atomic gas saturates. Dashed inclined lines represent 
“isochrones” of constant star formation efficiencies, indicating the depletion times = Ega^/^isFR of 10*, 10^ and 10*** years to consume all the 
gas. The dashed red isochrone marks a depletion time equal to the age of the universe, as one Hubble time (13.8 Gyr). A representative shift of the 
HaR-1 marker for Eh, is drawn, to show the "effective" molecular gas density at which stars would formed in this region if we consider a beam 
correction factor of +2.68 in log space, to convert our 22" beam to a ~1" beam as in the H„ data. 


in a +2.68 shift in log space, and is a representation of the real 
gas surface density at which stars would be formed, but always 
with the same r^ep. 

For Source-5 we have made the assumption that the HI and 
Ha emission are spatially correlated, i. e that they both belong to 
the gas plume associated to NGC4388. We know that this is true 
for HI, since we used the data from [Oosterloo & van Gorkom| 
( 2005| l and they probed the physical asociation of the HI gas 
plume to NGC4388. But this could not be the case for Ha, since 
we used the data from Kenney et al. (2008 1 and their Ha map 
show that this source could be associated to M86 instead, when 
Ha is considered. 

From Figure we can see that our sources have extremely 
low SFRs in comparison with the nearby spiral galaxies, and 
are only comparable with the most outer clumps in the H„/X- 
ray tail of the ISM stripped galaxy ES0137-001 in the Norma 
Cluster (Jachym et al. 2014|l and the XUV disk galaxies from 


[Watson et al.| ( |2014[ ) and Dessauges-Zavadsky et al.| ( |2014| l. We 


obtain depletion times significantly large. For example, HaR-1 
and HaR-2 have Tdep values of 2.2x10'^ and 1.2x10^' years re¬ 
spectively, to consume all the amount of gas present (HI+H 2 ). 
These values transform into 1.6x10'^ and 1.0x10^' years if we 
consider only the atomic gas component, and into 5.8x10" and 
1.6x10'° years if we consider only the molecular one. Such val¬ 
ues are quite large in comparison with the typical Tdep of ~2 Gyrs 
for spiral galaxies, and are even larger than a Hubble time by up 
to 2 orders of magnitude. In Table depletion times of H 2 that 
can be calculated are listed in column 7. 

The extremely low SF efficiency of our sources seems to fall 
off the linearity of the S-K relation for typical spiral galaxies at 
higher gas densities, a result previously reported in other low gas 
density environments, such as XUV disk galaxies ( [Dessauges^ 
[Zavadsky et al. |2014| l. [Watson et al. ( 2014| l presents a differ¬ 
ent conclusion for their results in XUVs, with a typical SFR in 
agreement with the S-K linear regime, but they take into account 
the 24//m emission in the SFR, and neglect the contribution of 
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heavy elements in the Zhj ■ We see that when we correct by these 
differences to make their data analytically compatible with ours 
(i.e neglect the 24jum emission and correct for heavy elements), 
their data points in the S-K plot are comparable to ours. 


4. Summary and Conclusions 

From our molecular cloud and star formation study in the tail 
north of NGC 4388 in Virgo, we can draw the following conclu¬ 
sions: 


1. CO(l-O) and CO(2-l) observations were done with the 
IRAM 30-m telescope in a total of 11 targets all along the 
ram-pressure stripped tail northeast NGC4388 in the Virgo 
Cluster, in order to probe the presence of molecular gas un¬ 
der extreme conditions. Such targets were selected for having 
strong peaks of HI and Ha emission. 

2. Four of such positions showed CO detections, and 3 of them 
concentrated in the HaR-2 region, at a distance of ~70 kpc 
of NGC4388, where molecular gas in unexpected. Given the 
large distances of these sources to NGC4388, it is not likely 
that the molecular gas was stripped from the galaxy, and 
must have formed in situ from the HI gas plume. 

3. Gaussian line profiles were fitted to the spectra of the detec¬ 
tions, finding a range of velocity dispersion between 12 and 
35 km/s. The CO(l-O) line profiles were used to estimate 
molecular gas masses and surface densities. The amount of 
molecular gas in these 3 regions (HaR-1, HaR-2 and HaR-2- 
4) is very low ( between 0.7 and 2.4 xIO^Mq), and their H 2 
surface densities between 0.2 and 0.9 Mopc^^. These values 
are well below the HI-H 2 threshold, where the gas is mainly 
atomic and very little is know about the SFR at such low gas 
densities, hence the importance of these detections. 


4. Using complementary data from|Yagi et al. (2013|l and Ken- 


|ney et al.l (120081) for Ha and from|Oosterloo & van Gorkom 
P005| l for HI, we computed Ssfr and Ehi to plot, in com¬ 
bination with Sh 2 , a S-K relation. Our sources show an ex¬ 
tremely low SFR (up to 2 order of magnitude lower than for 
typical spiral galaxies). For example, HaR-1 and HaR-2 have 
total gas depletion times of 2.2x10^^ and 1.2x10" years re¬ 
spectively. If we consider just the molecular gas component, 
these depletion times are 5.8x10" and 1.6xl0'° years. Fur¬ 
thermore, Source-1 and HaR-2-4 have H 2 depletion times 
even greater than 1.4 xlO" and 3.7 xlO" years respectively. 
These high values of depletion times exceed by far a Hubble 
time, thus indicating that this molecular gas will not form 
stars eventually, but remain in a gaseous phase and later join 


the ICM. 


5. From Figurej^we can see the linearity between the SFR and 
the gas surface density, so clear at high gas surface densities 
(> 9 Mopc^^) for normal spiral galaxies, cannot be extrap¬ 
olated to lower densities, below the HI-H 2 threshold, where 
the star formation is extremely inefficient, and the molecular 
gas, even though present, does not necessarily form stars . 
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